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We investigate the order-disorder transition line from a Bragg glass to an amorphous vortex 
glass in the H — T phase diagram of three-dimensional type-II superconductors taking into account 
both pinning-caused and thermal fluctuations of the vortex lattice. Our approach is based on the 
Lindemann criterion and on results of the collective pinning theory and generahzes previous work 
of other authors. It is shown that the shapes of the order-disorder transition line and the vortex 
lattice melting curve are determined only by the Ginzburg number, which characterizes thermal 
fluctuations, and by a parameter which describes the strength of the quenched disorder in the flux- 
line lattice. In the framework of this unified approach we obtain the H — T phase diagrams for both 
conventional and high-Tc superconductors. Several well-known experimental results concerning the 
fishtail effect and the phase diagram of high- Tc superconductors are naturally explained by assuming 
that a peak effect in the critical current density versus H signalizes the order-disorder transition 
line in superconductors with point defects. 
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I. INTRODUCTION 

In type-II superconductors one often observesEI a 
peak effect (or fishtail effect) in the critical current den- 
sity measured as a function of the applied magnetic field 
if at a fixed temperature T or as a function of T at 
fixed H. In conventional low-Tc materials this peak ef- 
fect mainly occurs at ^magnetic fields H near the upper 
critical field Hc2{T)Wb In the high-Tc YBaCuO crystals, 
the line of the maximum critical current density, Hp{T), 
frequently lies essentially below the irreversibility line in 
the H — T plane ,B~El and in sufficiently perfect crysljsJfi ji 
exhibits nonmonotonic behavior with temperature£3il3 
In these perfect crystals the line of maximum current 
density approaches the flux- line melting line, HrJJ^), 
approximately at the so-called upper critical poinlO at 
which the melting line terminates. When the oxygen 
deficiency 5 in YBa2Cu307_5 increases or the crystal 
becomes less perfect, the end point tends to the super- 
conducting transition temperature Tc at zero magnetic 
field, while_£/p(r) becomes a monotonically decreasing 
function.E2J't3 It is also important to note that at a fixed 
oxygen concentration the fishtail effect can disappear in 
pure YBaCuO crystals if the distribution o£-the oxygen 
vacancies over the sample becomes uniformJla 

At present the origin of the peak effect in low-Tc 
and high- Tc superconductors is commonly associated 
with thft-.proliferation of dislocations in the fluxJipe 
lattice.EJ"aE3"E3 At this first order phase transitionBlij"E3 
induced by quenched disorder in the vorte^ system, a 
transformation of a quasiordered Bragg glasscj into a dis- 
ordered amorphous vortex phase occurs. Although dif- 



ferent criteriaBliHil are used for determining the exact 
position of this transition on the peak-shaped dependence 
of the critical current density on H , they all lead to qual- 
itatively similar H — T phase diagrams, and, for definite- 
ness only, we shall imply below that the phase transition 
corresponds to the line of the maximum critical current 
density, Hp{T). 

A description of this order-disorder phase transition 
in high- Tc superconductors was proposed in Refs. p^~p6| 
using the Lindemann criterion. It was implied in these 
papers that the nature of the order-disorder phase tran- 
sition is different from the vortex lattice melting transi- 
tion, but at the critical point tpth phase transition lines 
merge. Recently it was refineifS that the upper critical 
point does not generally coincide with the point where 
the order-disorder line reaches the melting curve, and 
thus the melting line has a portion beyond the intersec- 
tion point. However, the following should be noted: The 
- 2^ for the disorder-induced transitiptt 



results of Refs. 24 



were obtained in the regime of single vortex pinningfZl 
when the Larkin pinning length Lc is less than Lq = ea 
where a = ($o/-ff)^^^ is the spacing between flux lines, 
$0 is the flux quantum, and e = Xab/K < 1 is the 
anisotropy of the superconductor {Xab, Ac are the Lon- 
don penetration depths in the ab plane and along the 
c axis, respectively). When T increases, the length Lc 
should exceed ea at some temperature which lies on the 
boundary of the single vortex pinning regime. At higher 
temperatures the disorder was completely neglected in 
Refs. |2^-p^, and only the melting line of the ideal lattice 
was derived. Thus, the behavior of the order-disorder 
line was not actually investigated in the high tempera- 
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ture region, and its connection with the melting hne was 
not established. Besides this, it has remained unclear 
why the proliferation of dislocations in the vortex lattice 
of high- Tc and low- Tc superconductors leads to different 
phase diagrams. 

In the present paper we consider the order-disorder 
transition line in the high temperature region and obtain 
the point where the vortex lattice melting and the order- 
disorder transition lines merge. It turns out that for a 
given model of the vortex pinning, the resulting H — T 
phase diagram is determined only by the Ginzburg num- 
ber Gi, which characterizes the thermal fluctuations, and 
by a parameteiEZi jc(0)/jo(0) that describes the strength 
of the quenched disorder in the flux-line lattice at T = 
(jo is the depairing current density and jc is the critical 
current density in the single vortex pinning regime; both 
are in the ab plane). For different values of these pa- 
rameters, phase diagrams are obtained which are similar 
to those observed in experiments for low- Tc and high- Tc 
superconductors. Thus, the results of this paper provide 
a unifled approach for analyzing H ~ T phase diagrams 
of various superconductors. 

In this paper we consider only magnetic flelds exceed- 
ing considerably the lower critical field Hd and thus dis- 
regard the reentrant behavior of the melting transition 
and do not distinguish between the magnetic field H and 
the magnetic induction B. Besides this, we deal only 
with anisotropic three-dimensional superconductors, ne- 
glecting completely the decoupling of the superconduct- 
ing layers. We also assume that H is directed along the c 
axis. This assumption simplifies the analysis of the prob- 
lem, though our final Eqs. (19) to (25) are valid for any 
direction of the magnetic field. 

II. LINDEMANN CRITERION 

We begin with simple estimates which show that the 
Lindemann criterion does define the condition for pro- 
liferation of dislocations in the fiux-line lattice at the 
order-disorder transition. Consider a dislocation net- 
work in the lattice. Let a unit cell of this network have 
the dimensions Rd and in the transverse and lon- 
gitudinal directions to H, respectively. A comparison 
of tilt and shear elastic energies yields that Ld/Rd ^ 
[c44(l/i?d, l/id)/c66]"^^^ > 1 where cge and C44(fcx,fc||) 
are the_shear and nonlocal tilt moduli of the flux-line 
lattice.Lj The energy cost for the creation of a disloca- 
tion cell is of the order of 

Ed ~ e^Ld 

where Eq — ('i'o/47rAaf,)^ and \ab is the London penetra- 
tion depth for currents in the ab plane. On the other 
hand, the elastic energy in the volume R^d^d is estimated 
as 

Eel ~ CGQLdU^{Rd, Ld) 



where (R, L) is the correlation function determining 
the relative displacement of points in the lattice with 
quenched disorder, 

u(i?,i)^([u(i?,i)-u(0,0)]2)i/^ 

Here u is the transverse displacement of a flux line, (...) 
means averaging over both thermal and quenched disor- 
der, the first coordinate R in u{R, L) indicates the posi- 
tion of the flux line in the plane normal to the applied 
magnetic field, while the second coordinate L defines the 
position of a point on the flux line. Comparing Ed and 
Eel with account of cge £o/a^i one arrives at the con- 
clusion that a dislocation network can exist in the lattice 
if u{Rd, Ld) > a. In other words, Rd should be greater 
than the so called positional correlation length Ra within 
which typical relative vortex displacements are of the or- 
der of the lattice spacing a. However, this is only a neces- 
sary but not sufficient condition for the existence of dis- 
locations. Displacements generated by the dislocations 
facilitate a better adjustment of the vortex lattice to the 
quenched disorder. The relative deformation of the lat- 
tice produced by the dislocations is of the order oia/Rd- 
Therefore, the smaller Rd is the greater is the gain SEpm 
in pinning energy Epin- Thus, the network first appears 
at the smallest possible Rd, and we arrive at the result 
Rd ~ Ra (and Ld ~ La) obtained in Ref. The rel- 
ative magnitude of the gain, S Epin/ Epin is determined 
by the ratio a/ Ra- Hence, for this magnitude to become 
of the order of unity, Ra/a should decrease to a certain 
constant C, 

a 

This criterion for the appearance of dislocations in the 
flux-line lattice was obtained in Refs. p2|^$| (see also 
Ref. IJ) , and is equivalentEJ to the conditioncJ: 

u^a,0)^cla^, (1) 

where is the phenomenological Lindemann constant. 
This immediately follows from the fact that the ratio of 
w(i?a,0) to u(a,0) [i.e., a/u(a, 0)] is a functiaji of Ra/a,- 
Finally, since u{a,0) = u(0, Lq) at Lq = eaj23 one more 
form of the Lindemann criterion exists: 

u\0,Lo)^cW- (2) 



It is just this form that was used in Refs. 




Strictly speaking, the values of the constants C and 
cl may depend on whether the order-disorder transition 
occurs in the single vortex pinning region or in the region 
of bundle pinning. However, to understand the essence 
of the matter, we shall use the simplest approximation: 
Cl will be considered as the same constant for the various 
regimes of pinning. 
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THE ORDER-DISORDER LINE. SIMPLIFIED 
APPROACH 



III. 



As well-known,EiHi!l thermal fluctuations of the flux- 
line lattice lead to a smoothing of the pinning potential 
and thereby affect the pinning. This thermal dcpinning is 
especially important for high-Tc superconductors. How- 
ever, to elucidate possible types of the order-disorder 
transition line, in this section we completely disregard 
the thermal fluctuations. The influence of the thermal 
depinning on the order-disorder line will be analyzed in 
Sec. IV. 



A. Region of single vortex pinning 

As has been mentioned above, the order-disorder line 
Hdis {T) was studiedEJ^Ej inside the single- vortex pinning 
regime where the Larkin pinning length Lc is less than 
Lq = ea. Siace Lq > L^., formulas of the random mani- 
fold regimaH3 for a single vortex are applicable to calcu- 
late the displacement correlation it(0, Lq), 



u(0,Lo) «e(io/ic) 



C 



(3) 



where ^ is the coherence length in the ab plane and C 
is the roughness exponent for a flux line. In Ref. p6| 
the value C ~ 3/5 was used, while C ~ 5/8 in Ref. |24 
Inserting Eq. (3) into Eq. (2), we obtain after simple ma- 
nipulations: 



Hdis — 



(4) 



where a = 2C/(1 — C) « 3. Eq. (4) coincides with the 
appropriate formulas of Refs. P4|-p6|. For Eq. (4) to be 
self-consistent, it is necessary to verify that Lq > Lc 
at H = Hdis or in other words Hdis < Hsv where 
Hsv — ^oc^/f-^c the boundary of the single vortex pin- 
ning regimeHj This condition yields: 



(5) 



If the inequality (5) is not fulflUed, Eq. (4) is not valid 
to describe Hdis- 

The parameter eS,/Lc generally depends on the tem- 
perature T. For example, according to simple estimates 
given in Appendix A, Eqs. (A3), (A7), it decreases with 
T and reaches zero at T = Tc- Moreover, its decrease 
becomes especially pronounced if the thermal depinning 
is taken into account. [In this case the single vortex col- 
lective pinning length Lc increases sharplyEZi when T ex- 
ceeds the characteristic pinning energy T^p{T), see Ap- 
pendix A.] Thus, even if €£_{0) / clLc{0) > 1, the order- 
disorder line Hdis {T) reaches the boundary of the single 
vortex pinning regime, Hs^{T), at some temperature Ti 
defined by the condition 



em) 

Lc{Ti) 



CL, 



(6) 



and at T > Ti, Eq. (4) fails. 

It is worth noting that the parameter tS,/ Lc appearing 
in Eqs. (4) - (6) and formulas given below characterizes 
the strength of the disorder in the flux-line latticeEZl and 
is expressed through the critical current density jc in the 
single vortex pinning regime. 



Lc 



1/2 



(7) 



where jo is the depairing current density. 



B. High temperature region 

At T > Ti the order-disorder transition line lies above 
Hsv {T) ■ In this region of the H — T plane ppall-bundle 
and large-bundle regimes of pinning occur. cZI Hence the 
transverse collective pii iruna I pngth Rc exceeds a, and to 
find M(a,0), the resultsE3SE3 may be used, which were 
obtained within the framework ofpthe perturbative ap- 
proach of Larkin and Ovchinnikov.E3 We have: 



u\a,Q)^e{LG/Lcf 



1 



1 



3/2 



(8) 



where Lq = ea, h = H/Hc2, hsv = Hsv/Hc2, the upper 
critical field Hc2 = ^o/'^tt^^, and Lc is the single vortex 
collective pinning length. Note that Eq. (8) differs from 
formula (4.17) of Ref. ^ by the last factor containing h 
and hsv This factor takes into account the possibility 
that 1 — /i is small; in Ref. ^ the correlation function 
u^{R, L) is given without taking account of this possibil- 
ity. The origin of this factor is the following. The quan- 



tity is proportional to nfl^^^\ab/Hc] 



1, 



one has \ab oc (1 

fpin ^0 



hy 



When h 

C66 oc (1 - /i)^ (see Ref. |2|||), 
while /pj„ cx £q cx (1 — h)"^ (see Ref. ^ and also Ap- 
pendix A). The combination of these factors gives Eq. (8), 
in which the additional constant factor (1 — hsv)^^^ has 
been introduced to provide a smooth crossover of this ex- 
pression to the appropriate formula for the single vortex 
pinning regime at h = hsv 

Inserting formula (8) into Eq. (1), we obtain an equa- 
tion for hdis = Hdis/Hc2, 



hdisi^ 



^dis } 



\clLc 



(1 



(9) 



where hsv — Hsv/Hc2 — 27r(e^/Lc)^ and the right hand 
side depends only on the temperature. Note that at 
T = Ti, when e£,/Lc = Cl, one has hdis{Ti) = 27rc| = 
hsv{Ti), in agreement with Eq. (4). A simple analysis 
shows that 27rc^ should be greater than 0.25 but less than 
1 (i.e., 0.2 < Cl < 0.4) for Eq. (9) to have a solution at 
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T > Ti. If CL < 0.2, the order-pdisorder line terminates 
at T = Ti, which is impossible .Ej On the other hand, if 
cl > 0.4, one finds from Eq. (9) that hdis > 1 — hsv, 
i.e., the root of the eqjiation Ues in the upper region of 
single vortex pinning23 where Eq. (9) is not vahd. For 
this reason, in the following we assume the conditions 
0.2 < Cl < 0.4 to be fulfilled and, for definiteness, take 
Cl — 0.25 in the subsequent calculations. 



C. Types of phase diagrams 

We begin the analysis of phase diagrams with the case 
D>CL 

where D = e^(0)/Lc(0) is the value of the parameter 
e£,/Lc at T — 0. Figure 1 shows the line Hdis{T) calcu- 
lated by solving Eq. (9) (for Hdis > Hsv) and Eq. (4) 
(for Hdis < Hsv)- In the construction of this figure, 
as well as in all examples below, we use Hc2(T) = 
Hc2i0)[l - (T/Tc)^] and a = 3 (i.e., C = 3/5). Besides 
this, taking into account the formula (A13) of Appendix 
A, we employ the following simple approximation for the 
parameter e^/Lc- 



Lc[T) 



Dgoit) 



with 



ffo(t) = (l-t^)^/^ 



(10) 



(11) 



where t = T/Tc- The increase of Hdis and its subsequent 
maximum are seen in the vicinity of Ti. As to hdis, this 
normalized quantity increases monotonically above Ti. 
When 1 — hdis ^ 1, an approximate solution of Eq. (9) 
is: 



HdisjT) 
Hc2{T) 



1/3 



(12) 



This formula shows that and how Hdis{T) approaches 
Hc2{T), which is also seen in Fig. 1. Interestingly, ac- 
cording to this formula, the order-disorder transition oc- 
curs outside the upper region of single vortex pinningtj, 
but its position correlates with the boundary of this re- 
gion: [1 - hd^s{t)]/[l - hsv{t)] = (27rc2)-2/3 > i. ft 
should be also noted that in the case under study (i.e. 
when the thermal fluctuations ara-negligible) the mean- 
field Hc2(T) practically coincidesEll with the melting line 
Hm{T). Thus, we obtain for D > cl that in the high tem- 
perature region a peak effect occurs near the melting line, 
while with decreasing T the position of the peak in jc{H) 
shifts downwards from this line. This situation iSpEcmj- 
niscent of that of perfect high- Tc superconductors.E3"llj 
In this context it is also useful to note the following: 
The density of the dislocations in the vortex liquid is es- 
sentially higher than in the disordered vortex solid phase 



near the order-disorder transitionE3 However, if in the 
H — T plane the order-disorder transition occurs suf- 
ficiently below the melting line of the clean supercon- 
ductor, then in the disordered solid phase, at the field 
corresponding to the melting transition, the density of 
the dislocations generated by the quenched disorder may 
become of the order of the density characteristic for the 
liquid phase. In this case the melting transition disap- 
pears. In other words, the melting line Hm{T) terminates 
when Hdis(t) deviates from it appreciably. 

If the strength of the disorder is sufficiently small, 

the order-disorder line lies entirely outside the region of 
single vortex pinning. Fig. 2, and is described by Eq. (9) 
at any T < Tc, while Eq. (4) is not valid at all. In 
the special situation when D is markedly less than cl, 
a peak effect occurs near Hc2{T) and its position in the 
H — T plane is approximately given by Eq. (12). In this 
case the resulting phct^e diagram looks like that of low- 
Tc superconductors.ErO The transition from one type of 
phase diagram to the other occurs when D = cl- 

It has been assumed in this section that the parame- 
ter e^/Lc decreases with increasing T. However, the ST^. 
pinning (due to spatial variations of T^) leads to an in- 
creasing function go{t), see Eq. (A14) in Appendix A. In 
this case, if D > cl, the formula (4) remains valid up to 
Tc- But if Z? < Cl, a temperature Tq exists, determined 
by the condition 



Lc{To) 



Cl 



and at T < To equation (9) should be used, while at 
T > Tq formula (4) holds. In other words, we have a 
situation which is opposite to that described above. In 
Fig. 3 the order-disorder line is shown for the case when 
go{t) is given by Eq. (A14). Note that in this case, ac- 
cording to formula (4), one has Hdis cx; (1 — i^)^/^ in the 
high temperature region of the phase diagram. Thi s se -, 
suit qualitatively agrees with the measurementsETBllilM 
on YBa2Cu307_5 crystals when 6 is not small, or when 
the crystals are not too perfect. 



IV. THE ORDER-DISORDER LINE WITH 
ACCOUNT OF THERMAL FLUCTUATIONS 

It is well-known that thermal fluctuations of the flux- 
line lattice play an important role in high- Tc supercon- 
ductors. In particular, for this reason the flux-line lattice 
melts essentially below the mean-field Hc2 line. In this 
section we study the influence of thermal fluctuations on 
the order-disorder line. 

Thermal fluctuations lead to a smoothing of the pin- 
ning potential and thus increase the Larkin length Lc- 
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The length Lc specifies the boundary of the single vor- 
tex pinning region and it enters the key parameter of the 
collective pinning theory, e^/Lc. Here we reserve the no- 
tation Lc for the true length renornialized by the fluctua- 
tions, while the Larkin length defined without account of 
the fiuctuations will be denoted below as L°. Note that 
just has been used in Sec. Ill, and just this quantity 
is described by Eq. (10). Apart from increasing Lc, the 
thermal fluctuations of the flux-line lattice also modify 
the correlation function (3) as follow£3: 



u2(0,Lo) «r2(Lo/Le)2«, 



(13) 



where = + w|., and ut is the magnitude of these 
fluctuations, which depends on the temperature and 
on the magnetic field. It is implied in Eq. (13) that 
Lo = ea > Lc- As to th£ correlation function (8), the 
collective pinning theorycZl gives 



0\3 



u'ia,0)»e{Lo/L'i) 



1 - h 



(14) 



for Lq < Lc, i.e. at /i > hgy — 27r(e^/Lc)^. 

Since u{a,0) = w(0, Lq), the correlation functions (13) 
and (14) must coincide at Lq = Lc (or equivalently at 
h ~ hgy). This condition yields 



(15) 



In fact, formula (15) is an equation for Lc since we have 
the relationship Hg^ — ^qe^ / L^. This equation enables 
us to find Lc{t), and thus Hsv{t), self-consistently. 

To proceed, we have to estimate the magnitude of ther- 
mal displacements of the lattice, ut, relative to its equi- 
librium position. In the case of the ideal vortex lattice, 
Ut was calculated in many papers; see, e.g., Refs. 
This magnitude, as well as the correlation functions 
and (14), depends on the elastic moduli of the lattice. 
However, in deriving Eqs. (8), (14) the contribution asso- 
ciated with the compression modulus cn was neglected. 
Hence, it is consistent to use the same approximation 
in the calculation of it^. This simplifies the appropriate 
formulacS for u'^, and we obtain 



Gi 



1/2 



(16) 



where h = H/Hc2{t), t = T/Tc, Hc2{t) = i?c2(0)(l - t^), 
Gi is the Ginzburg number. 



1 / Tc 
^'^2 Uc^(0)eV(0) 



which characterizes the strength of the thermal fluctua- 
tions, and He is the thermodynamic magnetic field of the 
superconductor. The function f(h) has the form: 



fih) 



2Pa [l + (l-f 5)^1^/^-1 
l-h c(H-c) 



(17) 



with c = 0.5[/3a(1 - ^)]^/^ and Pa = 1.16. 

The quenched disorder changes u'^. However, when 
the transverse collective pinning length Rc is consider- 
ably greater than a, the above result for the ideal lattice 
is a good approximation. This is due to the fact that the 
main contribution to ut results from the thermal fluc- 
tuations with short wavelengths {k± ~ 1/a), while the 
quenched disorder essentially distorts the lattice only on 
the scale Rc- Thus, we may use Eq. (16) in the case of 
the nonideal lattice if bundle pinning occurs. But in the 
single vortex pinning rcgime-thc influence of the disor- 
der is essential, and one ha£3 u^ (x Lc- To account for 
this result, we introduce an additional factor Lc/ea in 
the formula (16) and thus obtain 



Gi 

1 -i2 



1/2 



-1/2 



fih) 



(18) 



in the single vortex pinning region at ft. < hgy. 

Inserting expression (16) into formula (15) and using 
definition (10) for e£,/Lc, we obtain the following equa- 
tion for hsv = Hsy/Hc2- 



Gi 



1/2 



fiKvit))- (19) 



With increasing t the function hsv (t) reaches zero at some 
^dp < 1' ^^^d ^'i- (19) vahd at t < t^^. In the region 
t > the length Lc is infinite in size, and the single vor- 
tex pinning regime is absent, i.e., hsv(t) = a,t t > t 



dp' 



The value of is found by equating hgy to zero in for- 
mula (19): 



dp - r<a/2f(n\^°^^dp)l^ 



Gii/2/(o)- 



dp 



^2jl/2_ 



(20) 



It may be verified that the right hand side of Eq. (20) 
coincides (up to a numerical factor of the order of 
unity) with the dimensionless characteristic pinning en- 
ergy rjp/Tc, see Appendix A. Thus, in agreement with 
physical considerations,Ej we obtain that the essential 
renormalization of Lc occurs at such temperatures T that 

It should be noted that our result for Lc (Lc — > oo at 
t — > t^) differs in some respects from that presented in 
Ref. ^ where Lc increases exponentially at t t^^. How- 
ever, in the framework of our approximation, Hci = 0, 
we may consider Lc as infinite if it becomes of the order 
of A. Hence, the difference between the results is, in fact, 
small. But our approach provides the continuity of the 
correlation functions (13), (14) a.t h = hsv 

Equation (19) specifies the single vortex pinning (SVP) 
region existing at relatively low magnetic fields [0 < h < 
hsv{t)]. However, formulas (14), (16) enable one to find 
also the upper region in which the vortex system returns 
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to this type of pinning again. The appropriate equation 
results from the condition u(a, 0) = 
lowing form: 



and has the fol- 



(l-h) 



Gi 



1/2 



f{h) 



1 -i2^ 

27r {Dgo{t) f [1 - hs,,{t)], 



(21) 



As might be expected, at h — hgy this equation goes over 
into Eq. (19). However, in a certain temperature interval 
it has two additional real roots which form the boundary 
of the upper SVP region, (see Fig. 5 below). Here 
we do not consider this issue in detail, but qualitatively 
describe the effect of thermal fluctuations on the shape 
of the upper SVP region specified above without their 
account .E3 Although the softening of the vortex lattice 
near Hc2 is favorable for single vortex pinning, this soft- 
ness also leads to an increase of the thermal fluctuations 
ut, which reduces the strength of pinning. As a result, 
the upper region does not touch Hc2{t) except for the 
point t = (at t > and H = Hc2 we have ut = oo). 
Besides this, since ut increases with t, the upper region 
does not extend to Tc. Of course, one should realize that 
the part of the boundary of this SVP region lying above 
the vortex lattice melting curve has only formal meaning, 
since it does not account for the vanishing shear modulus 
at the melting. On the other hand, it is quite possible 
that the sharp melting transition disappears when the 
melting line enters the upper SVP region (or even before 
it; see Sec. V). 

Inserting expressions (13) and (18) in formula (2), we 
arrive at an equation for hdis = Hdis/Hc2 which general- 
izes formula (4): 



l + t 



Gi 



1/2 



f(hdis) 



1/(1-0 



n 2 \"/2 

2 2 / 27rc£ ^ 
^' hsvit) 



(22) 



This equation is valid in the single vortex pinning regime 
when hdis < hgv For example, if go{t) is a decreasing 
function of t, the order-disorder line lies in the single 
vortex pinning region and is described by Eq. (22) at 
t < ti. In the case D > the temperature ti is found 
from the equation: 



hsvih) = 27rc| 



l + < 



Gi 



1/2 



fihsvjt)) 

{hs.{t)Y'^ 



(23) 



which generalizes condition (6). At _D < cl the line is 
entirely outside this region, and thus one has ti = 0. 
If go{t) is an increasing function of temperature, more 
complicated situations can occur. 

Inserting formulas (14) and (16) into relation (1), we 
obtain the equation for hdis{t): 



hdisi^ — hdis)" 



l+t 



Gi 



1/2 



fihdi. 



i-ty {hd^sY'^ 

2^cl(^^^J[l-hsm', (24) 



which generalizes Eq. (9). This equation is valid in the 
bundle pinning region. 

Let us now present formulas for the melting line which 
is determined by the Lindemann criterion, = cj^a?, 
different from Eqs. (1), (2). This well-known empirical 
criterion based on the magnitude of the thermal fluctu- 
ations, was justified in Ref. ^ for the case of the ideal 
vortex lattice. According to this paper, different physi- 
cal mechanisms lead to the proliferation of dislocations at 
the vortex lattice melting and at the order-disorder tran- 
sition. While the disorder-induced transition is driven 
by an adjustment of the flux-line lattice to the disorder, 
the thermal melting is governed by the entropy gain as- 
sociated with the creation of dislocations. Hence, the 
Lindemann constant for the melting may, in princi- 
ple, differ from that used in Eqs. (22)-(24). However, 
for the sake of simplicity we take these constants cl as 
equal in the following analysis. Thus, if the melting line 
hm{t) = H„i/Hc2 does not intersect hsvit), it is described 
by the equation: 



Gi 



l-t^ 



1/2 



2ttcI. 



(25) 



But if the melting line lies below hgy (t) , an additional fac- 
tor {hm/hsvY^"^ should be inserted on the left hand side 
of Eq. (25), cf. Eqs. (16) and (18). Finally, we note that 
in expression (25), as well as in Eqs. (16), (18)-(24), the 
factor 1 — t^ represents the temperature dependence of 
the upper critical field, l — t^ = Hc2it)/ Hc2{0). Thus, if 
another form of this dependence is implied, the appropri- 
ate modification of all these equations is straightforward. 

Summing up, we may state the following: For given 
goit), equations (19)-(25) enable us to calculate the 
order-disorder line hdis{t), the melting line /im(i) and 
the boundaries of the single vortex pinning regime hsv (t) , 
h^P{t) with account of the thermal fluctuations. The 
function go{t) in Eq. (10) is determined by the pinning 
mechanism and by the temperature dependences of ^ and 
A, see Appendix A. 



V. ANALYSIS OF PHASE DIAGRAMS 

It is important to emphasize that for given goit) 
Eqs. (19)-(25) depend only on two parameters: the 
strength of quenched disorder, D = e^(0)/Lc(0), and the 
strength of the thermal fluctuations, Gi. In this sense, 
the equations and figures of Sec. HI correspond to the 
limiting case Gi — > 0. We consider now new features of 
the phase diagram which appear in the real situation of 
finite Gi. 
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A. Numerical results 



An example of the phase diagram in the case D > cl 
is shown in Fig. 4. Note that the order-disorder hne 
terminates at some temperature < 1. This termina- 
tion is associated with an increase of the thermal fluctu- 
ations and thus with an enhanced smoothing of the pin- 
ning potential when the order-disorder line approaches 
Hc2- InterestijiiglM, the end point lies near the so-called 
depinning lineEIlo where ~ S,"^ ■ Another new feature 
of the phase diagram is the intersection of the melting 
and the order-disorder lines. Thus, we obtain the point 
where both lines merge. It is seen from Fig. 4 that at 
this point the position of the critical current peak be- 
gins to shift sharply downward from thC|_mclting curve, 
in agreement with the experimental data£j1l3 Hence, as 
was mentioned in Sec. Ill C, the upper critical point for 
the melting is likely to occur somewhere nearby. The 
portion of the order-disorder line lying above the melt- 
ing curve has no physical meaning since the density of 
dislocations in the liquid plasise is already higher than in 
the disordered solid phase H3 Note also that the rise of 
Hdis{t) at T > Ti becomes considerably steeper than in 
the case Gi = 0. 

The intersection of the order-disorder line with the 
melting curve occurs also for D < cl- In this case, when 
Gi increases, the decrease of Hdis with T diminishes, and 
eventually Hdis becomes an increasing function of t, see 
Fig. 5. Therefore, if Gi ^ 10"^ (this value is typical 
for YBaCuO crystals), the temperature behavior of the 
order-disorder line is similar to that obtained for D > cl, 
and the phase diagrams in both these cases are of the 
same type. 

In Fig. 5 we also show the upper region of single vor- 
tex pinning. It is important to note that the intersection 
point of the order-disorder line and the melting line oc- 
curs clearly before the melting curve enters this region, 
and so at the intersection point the transverse collective 
pinning length Rc and the positional correlation length 
Ra (the dislocation spacing) both considerably exceed the 
flux-line spacing a. This fact justifies the assumed weak 
influence of quenched disorder on the vortex lattice melt- 
ing and the application of Eq. (24) up to the intersection 
point. 

If the parameter D increases, e.g., as a result of 
irradiation or of reduction of the oxygen content in 
YBa2Cu307_5 crystals (see Appendix A), then the 
order-disorder line shifts down, while the temperature of 
the intersection point (and thus the temperature of the 
upper critical point for the melting line) increases, Fig. 6. 
These r|es ulj . s . aij e in agreement with the experimental 
findingsta'Byiia for YBaaCugOy-s crystals. However, 
although the obtained results qualitatively describe the 
phase diagram of these crystals and its evolution with 
varying D, it should be emphasized that the intersec- 
tion point does not reach Tc at reasonable values of the 
disorder parameter D. 



Consider now the case when go{t) increases with t. 
This situation occurs in the model of 5Tc pinning, see Ap- 
pendix A. The appropriate phase diagrams are presented 
in Figs. 7-9. It is important to note the following: At suf- 
ficiently strong disorder, the line Hdis{t) monotonically 
decreases with temperature and practically reaches Tc- 
When the strength of the disorder decreases, the order- 
disorder line exhibits nonmonotonic behavior with t, and 
the temperature of the intersection point, ti, goes down. 
Thus, in contrast to the case of a decreasing function 
go{t), the model of STc pinning correctly reproduces all 
the featiirpsjof the experimental data for YBa2Cu307_5 
crystals.Hil3 In particular, the results presented in Fig. 8 
closely resemble the development of the order-disorder 
line with variations of the oxygen deficiency S, see, e.g.. 
Fig. 9 in Ref. |ll|. 

The results of this section may also provide an ex- 
planation of the findings obtained in Ref. 18. It was 



shown in this paper that the flux-line pinning in pure 
YBa2Cu307^5 crystals is mainly due to a nonuniform 
distribution of oxygen vacancies over the sample. Chang- 
ing the conditions of annealing of the sample, Erb et al. 
changed this distribution at a fixed 6. When the dis- 
tribution became more homogeneous, the fishtail effect 
disappeared. Since the spatial fluctuations in the density 
of oxygen vacancies strengthen the STc pinning, we get an 
additional confirmation of the hypothesis that this type 
of pinning plays the main role in not too perfect YBaCuO 
crystals, and so the fishtail effect exists up to high tem- 
peratures, as shown in Figs. 7, 8 for large D. When, 
after annealing, the crystal becomes more perfect, and 
the above-mentioned spatial fluctuations are reduced, the 
value of the parameter D appears to decrease consider- 
ably. (Moreover, it is conceivable that some other type 
of pinning begins to dominate.) Thus, we arrive at the 
situation when a fishtail effect is absent, at least in the 
region of not too high magnetic fields, as it follows from 
the data shown in Figs. 6 and 8 for small D. 

Finally, we briefly describe the case of small Ginzburg 
numbers which occurs for conventional superconductors. 
In this case the temperatures and ti are practically 
equal to unity, the melting line almost coincides with 
Hc2it), and the phase diagrams tend to those shown in 
Figs. 1-3. Thus, the results of the simplified approach of 
Sec. HI can be well applied to such superconductors. 



B. Approximate formulas. Discussion 

We now present analytical results which give some 
insight into the origin of the above-mentioned features 
of the phase diagrams. Let us start with the melting 
line. If [Gi/{1 - t2)]i/2 ^ 2ttcI, the normalized field 
hm = Hm/Hc2 is close to unity, and we obtain from 



Eq. (25) 



1 - K 



/i(l) 
27rci 



2/3 



+2/3 



Gi 



l-t^ 



1/3 



(26) 
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where 



hih)^{l-hf/'f{h) 



is defined through the function f{h) given by Eq. (17), 
with 1^(1) ~ 1-78. Equation (26) agrees with the strict 
resuliEll derived for clean superconductors without us- 
ing the Lindemann criterion. Formula (26) expressed in 
usual units {Hm = hmHc2) means that 

The right hand side of this expression is the width o, 
fluctuation region in not too small magnetic fields,! 
h ^ Gi. In the opposite limiting case when /i(0)[Gi/(l — 
^2-jji/2 27rc| [but the temperature is outside the criti- 
cal region for zero magnetic field, Gi < (1 — 
hm is small, and Eq. (25) yields: 



the field 



27rc| 
7^ 



Gi 



(27) 



wherfi 
resultl^^ 



« 2.34. In fact, this is the well-known 
for clean superconductors, 

if™ cx (1 - t)2. 



Note that for Eq. (27) to hold in a sufficiently wide tem- 
perature region, the Ginzburg number should not be too 
small. Finally, we point out an interesting feature of 
Eq. (25) that follows from our numerical analysis. The 
formula 



H„, = A{1 - t') 



2n7 



(28) 



with some constants A and 7 turns out to give a very 
good fit to Hm{t) determined by Eq. (25) in the wide 
temperature interval 0.5 < t < 0.98, see Fig. 10. The 
values of the parameters A and 7 depend on Gi, and the 
exponent 7 increases with increasing Gi. In particular, 
we find 1.24 < 7 < 1.59 when 0.001 < Gi < 0.01. In this 
connection it is worth noting that formulas of the type 
of Eq. (28) are widely used to approximate experimental 
data, and frequently an exponent 7 « 4/3 is found, which 
is characteristic for the fluctuation region of a 3D XY- 
type phase transition. We emphasize here that a good 
fit by such formulas does not necessarily mean the exis- 
tence of a large fluctuation region in zero magnetic fleld, 
but may result from the speciflc form of the expression 
[Eq. (25)] describing the melting line. 

In the case of STc pinning the melting line may lie in- 
side the region of single vortex pinning, i.e. the inequality 
hsv > may hold. This occurs if the parameter v (see 
Ref. P3 p. 1218), 



(2^)3/^1)3 



(29) 



is sufficiently large, v^^^ ^ 1. Although in this case the 
shape of the melting line requires a special investigation 



and is not discussed in detail here, it should be realized 
that the sharp melting transition may disappear in this 
region. 

Let us now analyze Eq. (24) which describes the line 
hdis{t) in the region h > hgv To understand the be- 
havior of this line, we note the following: The function 
fi{h) = (1 — hY^"^ f{h) decreases monotonically with in- 
creasing /i; its variation in the interval < /i < 1 is not 
large, /i(l)//i(0) ~ 0.76, and so if one considers fi{h) as 
constant, /i(/i) = /i(0) = /i « 2.34, this leads to a suf- 
ficiently accurate approximation in solving Eq. (24). In 
this approximation Eq. (24) becomes a quadratic equa- 
tion in the variable h^Jig{l — hdisY^'^ , and its solution has 
the form: 



"-die \^ 



Fn — Ft 



(30) 



where F^ and Ft are the following functions of temper- 
ature: 



2ci 



FT{t) = ht 



Gi 



1/2 



At t = we have Ft(0) = 0, and equation (30) goes over 
into Eq.(9). Thus, the effect of thermal ffuctuations is 
reduced to a renormalization of the right hand side of 
Eq. (9) [including the change of the function hsv{t) ac- 
cording to Eq. (19)]. Note that the magnitude of this 
renormalization (i.e., the ratio Ft/Fd) is mainly deter- 
mined by the parameter , see Eq. (29). Interestingly, 
in contrast to the case of the 5Tc pinning where v ap- 
pears as a result of the specific temperature dependence 
of go I Eq. (A14), we now conclude that the parameter v 
characterizes the relative strength of pinning of any type. 

It should be emphasized that the right hand side of 
Eq. (30) is real only if Fd > 2Ft. Thus, we find that the 
order-disorder line terminates at a temperature te which 
satisfies the equation 



FDite)^2FTite). 



(31) 



Taking into account Eqs. (16) and (31), it is easy to de- 
termine the ratio Ut/^^ at the end point of the order- 
disorder line, t = te, h = hdis{te)' 



FT{te)h'^'\l-h)-'/^ 



Frite 



Heiice, this end point lies on the so-called depinning 
lineE3 defined by the condition — S,^, see Figs. 4, 5. It 
can be also verified that t,, is always less than the tem- 
perature determined by Eq. (20). 

As has already been mentioned above, only the in- 
tersection point of the melting and the order-disorder 
lines has a physical meaning rather than the end point 
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of Hdisit). The temperature of this intersection point, 
ti, can be simply estimated using formulas (26) and (30). 
Then we arrive at the following equation for ti. 

AnclFoiU) = (1 + 2ncl)^FT{U), 

or explicitly, 

/i(l + 2^4)2 • ^^'^ 

This equation suggests that the temperature of the in- 
tersection point, ti, depends on the parameters Gi and 
D mainly through their combination Eq. (29). The 
data of Fig. 11 support this hypothesis, viz., at given 
go{t) the temperatures ti calculated for various Gi fall 
on the same curve. Of course, this prediction, as well 
as any quantitative conclusion based on Eqs. (19)-(25), 
requires an experimental verification since a number of 
simplifying assumptions were made above. In particular, 
cl was assumed to be the same constant for the melting 
and for the order-disorder transition in all the pinning 
regimes. On the other hand, according to Eq. (32), the 
dependence of ti on cl is relatively weak. 

Since ti specifies the width of the temperature region 
where the order-disorder transition exists, the data of 
Fig. 11 mean that this width is characterized by the pa- 
rameter ly. These data also shed light on the different 
behavior of Hdis{t) in Figs. 6 and 8 at large values of D. 
To elucidate the results of Fig. 11 and the conclusions 
made on their basis, let us analyze Eq. (32) in two lim- 
iting cases. If ^ 1, one obtains the following estimate 
from this equation: 

t, w 0.42(1 + 2TTciy^iy. 

In other words, we have a situation qualitatively simi- 
lar to that shown in Fig. 6 or in Fig. 8 for the smallest 
D. In this case the order-disorder transition occurs only 
at very low temperatures t < ti, while the portion of 
the melting line at t > ti has a large extension. Note 
that, according to the estimates of Gi and D presented 
in Ref. |[ such a situation must take place in pure crys- 
tals of 2H-NbSe2, which were investigated in numerous 
papers. This conclusion does not contradict the observa- 
tion of a peak effect in these crystals,B since a peak effect 
can signalize not only the order-disorder transition, as 
it has been— implied so far, but also the vortex lattice 
melting.C3 C3 Of course, the features of the peak effect 
may differ in these two cases. 

In the opposite limit, v ^ I, the order-disorder tran- 
sition reaches Tc, while the region of "pure" melting 
(ti < t < 1) is contracted. This limit just corresponds 
to conventional superconductors for which Gi 1, see 
Figs. 1-3. However, in the process of approaching this 
limit, the evolution of the phase diagram is different 
for the different models of pinning, i.e., in the cases of 
Eqs. (A13) and (A14). That is why Figs. 6 and 8 differ 



from each other at large values of D. If Eq. (A13) is 
valid, one finds from Eq. (32) that 

1-t^, ^ 1.5(H- 27rc|)z/-i/2. 

Thus, a large v is required for the temperature ti to reach 
unity, see Fig. 6. In the case of Eq. (A14) the right hand 
side of Eq. (32) already exceeds unity at a finite value ly 
{v ~ 1), and at such f the line Hdis{t) practically reaches 
the point T = Tc, H = 0, see Fig. 8 for the largest D. 



C. Upper critical point. 

Let us now discuss thcj-upper critical point on the vor- 
tex lattice melting curve.LlIn principle, two scenarios are 
possible. In the first oneE3 the upper critical point does 
not coincide with the intersection point of the melting 
and the order-disorder lines and is located at a higher 
magnetic field than this intersection point. In fact, in 
this case the intersection of the two different phase tran- 
sition lines occurs, and one of them (the order-disorder 
line) terminates at the intersection while the other one 
continues for some distance which seems to depend on the 
detailed position of the lines in the H—T plane. This sce- 
nario is based on the difference in the density of disloca- 
tions, p, generated at the melting and the order-disorder 
transition: at the melting one has p ^ a~^, while at the 
order-disorder transition p ~ R^^ a^^. The position 
of the upper critical point on the melting line is deter- 
mined by the condition that the density of dislocations in 
the disordered vortex solid phase reaches a value typical 
for the liquid (i.e., the dislocation spacing Ra calculated 
on the melting curve reduces to the vortex spacing a). It 
is this scenario that is assumed in our paper. 

In the second scenario the upper critical point coin- 
cides with the intersection point. In this case the above- 
mentioned point is ordinary rather than singular for the 
free energy of the vortex system, since a continuous phase 
transition separating the disordered vortex solid, phase 
from the vortex liquid does not seem to exist.CjC3 Thus, 
there is only one phase transition line which describes 
both the melting and the order-disorder transitions. This 
unified line originates from the melting curve of clean su- 
perconductors, gradually evolving from it as the strength 
of the quenched disorder increases. The "upper critical 
point" is then simply the point where on the melting line 
Ra reduces to a and the line experiences a bend. How- 
ever, for this scenario to occur, the calculated intersection 
point must lie in the single vortex pinning region. Other- 
wise, one obtains Ra > i?c 3> a at this point and returns 
to the first scenario. 

Our results argue in favor of the first scenario since in 
the framework of the used approximations the intersec- 
tion points never reach the upper region of single vortex 
pinning (see, e.g.. Fig. 5). Note that this scenario is also 
supported by recent numerical simulations.Eil 
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VI. CONCLUSIONS 

In this paper, using the Lindemann etiterion and re- 
sults of the collective pinning theorypO we have de- 
rived equations (19)-(25) which enable one to calculate 
the order-disorder transition line Hj^is(t) with account 
of both pinning-caused and thermal fluctuations in the 
whole temperature interval of its existence. The bound- 
aries of the single vortex pinning region, Hsv{t), Hl'-P{t), 
and the melting line Hm{t) are also found from these 
equations. The equations turn out to depend only on the 
Ginzburg number Gi which characterizes thermal fluctu- 
ations, on the strength of the quenched disorder, Z?, 



D = 



Lc{0) Vjo(O) 



1/2 



and on the function go{t) defined by the pinning mech- 
anism and by the temperature dependences of ^ and A, 
Appendix A. For example, the pinning mechanism con- 
sidered in Ref. ^ leacLsjo a go{t) described by Eq. (A13), 
while the STc pinning23 results in Eq. (A14). Using our 
equations one can analyze the phase diagrams of vari- 
ous superconductors. Moreover, our analysis in princi- 
ple allows to obtain information on the true form of the 
function go{t), i.e., on the mechanism of pinning in a 
superconductor. 

At small Gi we obtain phase diagrams typical for low- 
Tc superconductors. In this case the obtained results 
practically coincide with the results of the simplified ap- 
proach presented in Sec. III. We also analyze phase dia- 
grams with Gi ~ 10^'^ since such values of Gi are char- 
acteristic for the YBa2Cu307_5 superconductors. We 
consider both a typical case where go{t) decreases with 
t [Eq. (A13)] and a case where go{t) increases with t 
[Eq. (A14)] . In both these cases the obtained results qual- 
itatively describe the phase diagram of YBaCuO crystals 
and its evolution with varying Z?, which can be caused by 
irradiation or by reduction of the oxygen content in these 
crystals. However, we find that the model of 6Tc pinning 
[Eq. (A14)] more correctly reproduces all the features of 
the experimental data for YBaaCusOy-^ crystals. 

Our results support the ideacHl that the upper critical 
point on the melting curve of high- Tc superconductors 
does not generally coincide with the intersection point of 
the melting and the order-disorder lines (although the 
distance between them is possibly small). We predict 
the position of this intersection point in the iJ — T plane. 
We suppose that the temperature corresponding to this 
point is determined mainly by the parameter i/, Eq. (29), 
i.e., by a combination of the parameters D and Gi. 
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APPENDIX A: ESTIMATES OF THE 
PARAMETERS 



The results Qt|tiiis paper are based on the collective 
pinning theory.EIlE^ It is assumed in this theory that 
the disorder in the flux-line lattice is generated by point 
defects of size not exceeding the coherence length ^(0). 
Here, in the framework of three widely used models of 
pinning, we express the quantities e^/Lc and TJp through 
the characteristics of the point defects, viz., their concen- 
tration n and mean radius vq. These expressions enable 
us to get some idea of the temperature dependences of 
these quantities and to understand the changes of the 
phase diagram when_a. and vq are varied. 

In the first model,t3 pinning is due to the gain in the 
condensation energy when a vortex core is located at a 
defect. This gain €pin is of the order (iZ^/87r)47rrQ/3, 
where He is the thermodynamic magnetic field and vq 
the radius of the defect. Then one obtains the following 
estimates for the mean elementary pinning force exerted 
by one point defect /p™ ~ epi„/^: 



fpin 



sol J 



(Al) 



and for the single vortex collective pinning lengthE Lc = 

n-'/H^^eo/fp,nf^^ [with e = Xat/Xc, £o = ($o/47rAa6)2]: 



£4/3 



(A2) 



Using these estimates, we find the key parameter (7), 



Lr. 



3\l/3 



(A3) 



and the characteristic pinning energycZI TJ^ described by 
the expression ^(/^^^nLc^^)!/^, 



dp 



(A4) 



A different mechanism of pinning is due to the scat- 
tering of quasiparticles by the defect as calculated in 
Ref. p2l A scattering center facilitates deformations of 
the order parameter up to distances of the order of the 
zero-temperature coherence length ^(0). Hence, it is en- 
ergetically advantageous for a vortex core to sit at a scat- 
tering center. This mechanism leads to the pinning force 
fpin ~ (ffc/87r)ro^(0)/^, and we obtain 
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fpin 

Lc 



roxm 



[roj 



2/3 



T, 



dp 



(A5) 

(A6) 

(A7) 
(A8) 



Strictly speaking, formulas (A5)-(A8) are valid if 
7rrQ^(0)B-jC 1. This type of pinning is sometimes called SI 
pinning .E3 Note that equations (A5)-(A8) and (A1)-(A4) 
have the same dependences on n, e and temperature. 
They thus describe two different contributions to one pin- 
ning mechanism, with the contribution of Eqs. (A5)-(A8) 
always dominating when ro < £.(0)- 

If the pinning centers in YBa2Cu307_5 are clusters of 
oxygen vacancies, it is useful to keep in mind that 



3 r 

—rnU = CO, 



(A9) 



where the constant c equals 2 if these clusters are formed 
by YBa2Cu306.5 phase. Thus, when ro decreases at a 
fixed (5, the parameters e^/Lc and TJ^ also decrease. 

In high- Tc superconductors the pinning can be due to 
spatial fluctuations in the density of the oxygen vacan- 
cies, which results in variations of Tc over the sample. 
One has the following estimates for this STc pinningiEZI 



1/3 



(AlO) 
(All) 



where n is the density of oxygen vacancies. This pinning 
can occur when ne^(O)^ > 1. 

The above equations enable us to estimate the temper- 
ature dependence of the quantities e^/Lc and TJ^. For 
this purpose, we insert the expressions 



m m 



T2 



-1/2 



(A12) 



in the appropriate formulas. In particular, it follows from 
Eqs. (10), (A3), (A7) that for SI pinning 



(A13) 



where t = T/Tc, D = e^(0)/Lc(0), while according to 
Eq. (AlO), one finds for STc pinning 



go{t) = {l^t^)-''^. 



(A14) 
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FIG. 1. The order-disorder line Hdis{t) (solid line) calcu- 
lated from Eqs. (4), (9), (A13) {51 pinning) for cl = 0.25 
and D/cl = 1.3. Simplified approach, without account 
of thermal fluctuations. The boundary of the single vor- 
tex pirming regime, Hav{t), is given by the dashed line, 
and the dotted line shows the mean-field upper critical field 
H,2(t) = H,2(Q)(l-t'). 
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FIG. 3. The hnes Hdu{t) and Hsv{t) as in Figs. 1, 2 but for 
5Tc pinning, Eq. (A14), for D/cl = 1.3 and 0.9. Now Hdis{t) 
and Hsv{t) cross at t = to for D/cl < 1- 
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FIG. 2. As Fig. 1 but for D/cl = 0.9 and 0.4. Note that for 
D/cl < 1 the order-disorder line Hdia{t) does not intersect 
the single vortex pinning boundary Hav{t). 
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calculated from Eqs. (22), (24), (A13) (5/ pinning), which ac- 
count for thermal fluctuations, for cl = 0.25, D/cl = 1.3, and 
for three values of the Ginzburg number Gi = 0.01, 0.001, and 
0.0001 (from top to bottom). Also shown are the boundary 
of the single vortex pinning regime Hsv{t), Eq. (19) (dashed 
line), the mean- field Hc2it) = Hc2(0)(l - t'^) (dotted line), 
the vortex- lattice melting line Hm{t), Eq. (25) (dashed-dotted 
line), and the depinning line Hdp(t), Eqs. (16)-(18) (thin solid 
line), along which one formally has = and where Hdis{t) 
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FIG. 5. As Fig. 4 but for different D/cl = 0.9. The dashed 

lines show the boundaries of both the lower (Hsv) and upper 
(Hgjj) regions of single vortex pinning, see Eq. (21) and text. 
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FIG. 6. The order-disorder lines Hdis(t) (bold solid lines) FIG. 8. The order-disorder lines Hdis{t) (solid lines) as in 

according to Eqs. (22), (24), (A13) (SI pinning) for cl = 0.25, Fig- (7) with (A14), but for many values D/cl = 0.6 to 2. 
Gi = 0.01 and several values D/cl = 0.4 to 1.6. The 
dash-dotted line is the vortex-lattice melting line Hm{t), 
Eq. (25), and the dotted line indicates i?c2(t) = i/c2(0)(l-t^). 




in Fig. 6 but for STc pinning, (A14), for D/cl = 1.3, 
1.1, and 0.9. The dashed lines are the single-vortex pin- 
ning boundaries Hsv{t), Eq. (19), the dash-dotted line is the 
melting line Hm{t), Eq. (25), and the dotted line indicates 

Ha2it)=Hc2iO){l-t''). 
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FIG. 10. The vortex- lattice melting-line Hm(t) from 
Eq. (25) (solid line) and its approximations Eq. (26) (crosses) 
and Eq. (28) (circles) with A = 0.82 and 7 = 1.31. The dashed 
line indicates Hc2{t) oc l-t^. Here cl = 0.25 and Gi = 0.002. 
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FIG. 11. The temperature U where the melting line and 
the order-disorder line cross, plotted versus the parameter 

u = {2tt)^^^ D^/Gi^^^. This temperature is found by calculat- 
ing the phase diagrams from Eqs. (19)-(25) with cl = 0.25 
and Eq. (A13) (right 3 curves) or Eq. (A14) (left 2 curves) for 
Gi = 0.01 (crosses), Gi = 0.001 (triangles), and Gi = 0.0001 
(circles). The interval of D/cl actually used for each curve 
can be calculated from Gi and the v interval. 
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